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ABSTRACT 

We present radio observations of GRB 980703 at 1 .43, 4.86, and 8.46 GHz for the period of 350 to 1000 days after 
the burst. These radio data clearly indicate that there is a persistent source at the position of GRB 980703 with a 
flux density of approximately 70 /iJy at 1.43 GHz, and a spectral index, (3 « 0.32, where F v oc v~P . We show that 
emission from the afterglow of GRB 980703 is expected to be one to two orders of magnitude fainter, and therefore 
cannot account for these observations. We interpret this persistent emission as coming from the host galaxy — 
the first example of a 7-ray burst (GRB) host detection at radio wavelengths. We show that emission from an 
AGN is unlikely, and find that it can be explained as a result of a star-formation rate (SFR) of massive stars (M > 5 
M Q ) of w 90 M Q /yr, which gives a total SFR (0.1M Q <M< 100M Q ) of w 500 M Q /yr. Using the correlation 
between the radio and far-IR (FIR) luminosities of star-forming galaxies, we find that the host of GRB 980703 is 
at the faint end of the class of Ultra Luminous Infrared Galaxies (ULIRGs), with LfiR ~ few x 10 12 L Q . From 
the radio measurements of the offset between the burst and the host, and the size of the host, we conclude that 
GRB 980703 occurred near the center of the galaxy in a region of maximum star formation. A comparison of 
the properties of this galaxy with radio and optical surveys at a similar redshift (z ~ 1) reveals that the host of 
GRB 980703 is an average star-forming galaxy. This result has significant implications for the potential use of 
a GRB-selected galaxy sample for the study of galaxies and the IGM at high redshifts, especially using radio 
observations, which are insensitive to extinction by dust and provide an unbiased estimate of the SFR through the 
well-known radio-FIR correlation. 

Subject headings: gamma rays:bursts - radio continuum:general - cosmology:observations - galaxies:starburst - 
stars:formation 



1. INTRODUCTION 

Recent studies of the properties and host galaxies of 7-ray 
bursts (GRBs) reveal some indirect evidence for the link be- 
tween GRBs and star formation. Optical measurements of the 
offset distribution of GRBs from their host centers appears to be 
consistent with the distribution of collapsars in an exponential 
disk, but inconsistent w ith the expected offset distribution o f 
delayed binary mergers (Bloom, Kulkarni & Djorgovski 2001 ). 
GRB 990705 is an illustrative example of this result since HST 
images revealed that the burst was situated in a spiral arm, just 
north of an apparent star forming reg ion (Holland et al. 2000; 



Vreeswijk et al. 1999) offers an excellent opportunity for 



Bloom, Kulkarni & Djorgovski 2001 ). The absence of optical 
afterglows from the so-called "dark GRBs" (Djorgovski et al. 
2001) points to the association of GRBs with heavily obscured, 
and possibly star-forming, regions. In addition, Galama & Wi- 
jers (2000) claim high column densities toward several GRBs 
from X-ray observations of afterglows 

Consequently, one of the pressing questions in the study of 
GRB host galaxies is whether they are representative of star- 
forming galaxies at a similar redshift. If they are, then the 
dust-penetrating power of GRBs and their broadband afterglow 
emission offer a number of unique diagnostics of their host 
galaxies: the obscured star formation fraction, the ISM within 
the disk, the local environment of the burst itself, and global 
and line-of-sight extinction, to name a few. 

GRB 980703, which has one of the bri ghtest (apparent mag 
nitude) hosts to date (R w 22.6 mag; Bloom et al. 1998 



detailed studies. The afterglow optical and near-IR (NIR) 
lightcurves exhibited pronounced flattening about 6 days after 
t he burst and this was attributed to an under lying bright host 
( [Bloom etal. 1998t |Castro-Tirado et al. 1999|; V reeswijk et al. 
1999). Djorgovski et al. (1998) undertook spectroscopic obser- 
vations of the host and obtained a redshift of 0.966. Using three 
different estimators ([Oil], Ha and 2800A UV continuum) of 
the star-formation rate (SFR), Djorgovski et al. (1998) inferred 
extinction-corrected SFR of 10 to 30 MQ/yr. 

Here we report radio observations of GRB 980703 covering 
the period 350-1000 days after the burst at three frequencies: 
1.43, 4.86, and 8.46 GHz. This burst has the distinction of be- 
ing followed up for 100 days; the previous record-holde r was 
GRB 970508 (445 days; |Frail, Waxman & Kulkarni 2000| ). The 
organization of the paper is as follows. We summarize the radio 
observations and data reduction in §|[ In §|| we show that the 
late time radio observations require a steady component over 
and above the decaying afterglow observations. We argue that 
this component is unlikely to arise from an AGN but is instead 
due to star-formation. In §|], we infer the SFR from the radio 
observations and compare and contrast this estimate to those 
derived from optical observations. Thanks to the high angu- 
lar resolution and accurate astrometry of radio observations we 
are able to derive an accurate offset between the burst and the 
centroid of the host, as well as constrain the size of the radio 
emitting region (§q). 
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2. RADIO OBSERVATIONS 

Very Large Array (VLA 3 ) observations of GRB 980703 were 
initiated on 1998, July 4.40 UT at 4.86 GHz. All observations 
were obtained in the standard continuum mode with 2 x 50 
MHz contiguous bands. We used the extra-galactic sources 
J2330+110, J0010+109 and J0022+061 for phase calibration 
and 3C48 (J0137+331) and 3C147 (J0542+498) for flux cali- 
bration. We used the Astronomical Image Processing System 
(AIPS) for data reduction. The fluxes presented here are the 
values measured on the images at the position of the source. 

Late-time observations (time after the burst, t ;> 350 days) 
were co-added over a period of a few to thirty days in order to 
increase the overall sensitivity of each detection. This is ap- 
propriate since the expected change in the flux density from the 
afterglow over a few days, several hundred days after the burst, 
is negligible relative to the associated errors in the measure- 
ments. A log of the late-time observations and the flux density 
measurements are summarized in Table |J and the lightcurves 
are shown in Figure |l| A summary of the early radio data, as 
well as broadband modeling is given in Berger et al. (2001). 



3. EVIDENCE FOR HOST GALAXY EMISSION IN THE RADIO 
REGIME 

From Figure [l], we see that the late-time (t > 350 days) radio 
lightcurves do not exhibit the customary power-law decay ex- 
pected of afterglows but instead show flattening. From the early 
broadband data we know that the afterglow spectrum peaked at 
freq uency, v m ~4x 10 12 Hz at t = 1.2 days (V reeswijk et al. 
1999). If the explosion was spherical then we expec t v m oc f" 3 / 2 
([Sari, Piran & Narayan 1998|; phevalier & Li 2000|). Thus the 



radio afterglow is expected to decay for t > 70 days after the 
burst. If the ejecta were collimated (opening angle, 6j), then 
we expect a more rapid decay, v m oc f~ 2 , once the bulk Lorentz 
factor, T, of th e flow falls below 6>,, T(t) < 9~ : l ( Sari, Piran 
& Halpern 1999). In this case, we expect the radio afterglow 
to start decaying at even earlier times, and the flux will decay 
faster relative to a spherical explosion. In either case, we expect 
the radio afterglow to decay by at least a factor of three over the 
time span under consideration, 350 < t < 1000 days. 

We can clearly see from Figure [j] that this decay is not taking 
place, and the flux instead remains constant over a period of ap- 
proximately 650 days. This behavior is similar to the flattening 
observed in the optical/NIR lightcurves of several GRBs (in- 
cluding GRB 980703), when the emission from the afterglow 
decays below the level of emission from the host galaxy. 

Furthermore, the afterglow spectrum is expected to be a 
power law, F v oc v~°, where (3 = (p— l)/2 and p is the power 
law index of the Lorentz factor distribution of the shocked elec- 
trons , N(^i)d~f oc r y~ p d^f for 7 > 7 m j n (S ari, Piran & Narayan 
1998). From the observations of many afterglows, we note that 
p is in the range 2.2-2.6 and thus we expect f3 ~ 0.7. How- 
ever, the observed spectral index in the range 1.43-8.46 GHz is 
much lower, (3 = 0.32 ±0.12. We thus conclude that there exists 
a steady source of emission other than the afterglow. 

One possible explanation for this component is emission 
from an active galactic nucleus (AGN). It has been noted in sur- 
veys of the Hubble Deep Field (HDF), its flanking fields, and 
the Small Selected Area 13 (SSA13) that approximately 20% 

3 The VLA is operated by the National Radio Astronomy Observatory (NRAO), 
by Associated Universities, Inc. 



of t he radio sources are AGN with spectra l indices of about 
0.3 ( |Richards et al. 1999|; [Richards 2000[ Barger, Cowie & 
Richards 2000). Windhorst et al. (1993) found a similar re- 
sult in their survey of two 7' x 7' fields at 8.44 GHz. Thus, 
there is a modest probability that the emission from the host of 
GRB 980703 is due to an AGN. 

We consider the AGN hypothesis unlikely based on the ra- 
dio data and optical spectroscopy. First, optical spectra of 
the source obtained by Djorgovski et al. (1998) show no 
evidence for an unobscured AGN: high-ionization lines such 
as Mg II A2799, [NeV]A3346, and [NeV]A3426 are absent, 
and the [OIII]A4959 to H/3 ratio is approximately 0.4, much 
lower than [OIH1/H3 > 1.3 for AGN ( Rola, Terlevich & Ter- 
levich 1997). Another way to discriminate between AGN and 
star-forming galaxies is to co rrelate the [OH] equival ent width 
(EW) with continuum color (Dressier & Gunn 1982). Kenni- 



cutt (1992) showed that AGN have redder colors for similar 
[Oil] EW, relative to normal galaxies. Using the spectrum pre- 
sented in Djorgovski et al. (1998) we evaluat e the color index , 
(41-50) = 2.51og[/„(5000A)/%(4100A)] ( [Kennicutt 1992[ ), 
and find it to be ± 1; an AGN with the same [Oil] EW would 
have a value ;> 0.3 (Kennicutt 1992). Finally, Rola, Terlevich 
& Terlevich (1997) found that for a sample of emission-line 
galaxies at z ~ 0.8, the color index between the continuum un- 
derlying the H/3 and [OII]A3727 lines is > 0.4 for all AGN in 
their sample. Using the spectrum of Djorgovski et al. (1998) 
we find that this color index is approximately zero. 

A second, but less persuasive argument against an AGN ori- 
gin is the apparent absence of significant radio variability over 
the 650 day monitoring period (see Figure ||). The radio cores 
of most, but not all, low-lu minosity AGN show variability ex- 
ceeding the observed levels ( ^alcke et al. 2000| ). 

We thus conclude that the radio emission seen from the host 
of GRB 980703 is unlikely to be due to AGN activity. How- 
ever, star-forming galaxies exhibit radio emission arising from 
their supernova remnants (SNRs) and HII regions. In the next 
section we show how the observations, the radio spectral index, 
and the optical spectrum, are consistent with the hypothesis that 
the radio emission is related to star formation. 

4. THE STAR-FORMATION RATE IN THE HOST GALAXY OF 
GRB 980703 

Star formation is traced by optical, far-IR, sub-mm, and radio 
emission. In the following we will use the radio data to estimate 
the SFR in the host galaxy of GRB 980703, and then compare 
the results with the SFR derived from optical indicators, and 
with radio surveys at a similar redshift range in order to place 
the host of GRB 980703 in a larger context. 

4.1. Star Formation Rate from the Radio Observations 

Using all the measurements in Table [j], we find the fol- 
lowing weighted-average flux densities for the host galaxy of 
GRB 980703: F uAA(l = 39.3 ± 4.9 ply, F vA%6 = 42.1 ± 8.6 pJy, 
and /v. 1.43 =68.0±6.6^Jy. From t he redshift of GRB 980703, 
z = 0.966 ( Djorgovski et al. 1998 ), and the cosmological pa- 
rameters Slo = 0.3, Ao = 0.7 and Hq = 65 km/sec/Mpc, we find 
that the luminosity distance to the burst is d^ = flU(l +z) 2 ~ 
2.1 x 10 28 cm, and the observed luminosity at each frequency 
is given by L„ = \itd\F v . The emitted luminosity is given by 
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L em y = L bs.v(v jv'YiX + z) , where (3 is the spectral index of 
the radio emission, v' is the rest frequency, and v is the observ- 
ing frequency. From the flux densities we find that the radio 
spectral index is (3 = 0.32 ±0.12, and thus the emitted lumi- 
nosity in each frequency is approximately 25% higher than the 
observed luminosity at the same frequency. 

Condon (1992) showed that the total luminosity is a combi- 
nation of synchrotron and thermal emission components, both 
directly related to the formation rate of massive stars via a 
simple relationship. Moreover, since the lifetime of massive 
stars is of the order of 10 7 years, and the lifetime of the syn- 
chro tron emitting electrons is of the order of 10 8 years (C ondon 
1992), the radio emission is an excellent probe of the instan- 
taneous SFR. Using the emitted luminosity at v 1 = 1 .43 GHz, 
L em (1.43) = (4.7 ± 0.6) x 10 30 erg sec" 1 , we find that the SFR 
of massive stars in the host of GRB 980703 is 



SFR(M>5M ); 



L m (1.43) 



--«9OM /yr. 

5.3 x 10 2 V^ z + 5.5x 10 2 VF 01 



GHz 



(1) 

Since both the thermal and non-thermal components are pro- 
portional only to the formation rate of high-mass stars, equa- 
tion [l] has to be modified by a factor which accounts for the con- 
tribution from stars in the mass range 0. 1-5 M . For a Salpeter 
IMF this factor evaluates to 5.5. We use the Salpeter IMF since 
it is already implicitly used in equation [l] for the mass range 
5-100 M . Thus, within this framework the total SFR is sa 500 
M Q /yr. 

We now compare the SFR and other characteristics of the 
host of GRB 980703 with those of other galaxies at a similar 
redshift. A radio survey of the Hubble Deep Field (HDF) and 
its flanking fields showed that the mean spectral index of the 
8.46 GHz selected sample is (/3x M ) = 0.35 ± 0.07 (R ichards 
2000). In a survey of two 7' x 7' fields with the VLA at 8.44 
GHz Windhorst et al. (1993) found for sources with a flux 
density <J 100 /iJy (hereafter, /iJy sources) a median spectral 
index, (3 med w 0.35 ±0.15, and Fomalont et al. (1991) found 
Aned ~ 0.38 for /iJy sources selected at 4.9 GHz. Thus, the 
host galaxy of GRB 980703 appears to be a normal /iJy source 
compared to sour ces selected at 4.9 or 8.5 GHz. In addition, 
it has been noted (Richards 2000) that the spectral index of ra- 
dio sources selected at frequencies larger than 5 GHz flattens 
from a value of approximately 0.7 for the mJy (F v ;> 1 mJy) 
population to 0.3 for /iJy sources. 

The reason for the flattening of the spectral index is a varying 
ratio of thermal bremsstrahlung to synchrotron emission. Su- 
pernova remnant shock acceleration of electrons results in syn- 
chrotron emiss ion, with a characteristic spectral index of m 0.8 
flCondon 199^ ). On the other hand, thermal bremsstrahlung 
emission from HII regions, excited by star formation, has a 
much natter spectral index, « 0.1. Thus, as the direct con- 
tribution from massive stars increases the spectra are expected 
to flatten from a value of 0.8 to 0.1. This is exactly the effect 
that is observed in the aforementioned surveys. 

Within the HDF and SSA13 Richards et al. (1999) identified 
radio sources with fluxes in the range 10-100 /iJy with bright 
disk galaxies with / « 22 mag. The I— K color for these galax- 
ies is approximately 2.5 mag. Bloom et al. (1998) find 21.9 
mag and I-K w 2.1 mag for the host of GRB 980703. Thus, we 
see from both the radio spectrum of the source, and the optical / 
mag and I-K color that the host galaxy of GRB 980703 has the 



characteristics of a typical star-forming radio galaxy selected at 
8.5 GHz. 

Figure || shows the total SFR in the host galaxy of 
GRB 980703 as compare d to sources in the HD F, its flanking 
fields, SSA13, and V15 ( fjaarsma et al. 2000| ) in the redshift 
range 0.85- 1.15. These fields have been observed to /iJy sen- 
sitivities at cm wavelengths, and the detected radio sources have 
been identified with optical sources for which the redshift was 
determined. We use the flux and spectral index measurements 
along with equation [j] and the correction factor to calculate the 
total SFR. It is clearly seen from the figure that the host of 
GRB 980703 is an average galaxy at z ~ 1 based on star for- 
mation. This conclusion meshes well with the comparison of 
the radio spectral index, optical / mag, and optical I-K color 
of the host of GRB 980703 to the same sample. 

4.2. Star Formation Rate from Optical and Sub-mm Data 

Djorgovski et al. (1998) used Ha and the 2800A UV con- 
tinuum to calculate a SFR of approximately 10 M /yr in the 
host of GRB 980703, after correcting for rest-frame extinction, 
Ay 0.3 mag. Sokolov et al. (2001) found a similar intrinsic 
extinction, Ay w 0.3-0.65, and based on template spectral en- 
ergy distributions found that the best model for the broadband 
optical spectrum is given by exponentially decreasing star for- 
mation with an extinction-corrected SFR of 20 M /yr. 

Clearly, the SFR derived from optical indicators is much 
lower than the value from radio measurements, even after cor- 
recting for extinction. This result is part of a general trend that 
has been obser ved in galaxies with SFR > 0.1 M^/yr (H op- 
kins et al. 2001). Hopkins et al. (2001) propose dust redden- 
ing dependent on SFR as the solution to this problem, and we 
therefore expect a much better result if we use their prescrip- 
tion. Extending their correlation to SFR1.43 w 500 M /yr, we 
find that the predicted observed SFR from Ha is approximately 
70 M /yr. This value is still much higher than the measured 
SFR. In fact, in the Hopkins et al. (2001) sample the optically- 
derived SFR rarely exceeds 10 M /yr, while the radio-derived 
values go up to several hundred M /yr, indicating that the op- 
tical emission does not trace the entire star-forming volume. 
Thus, the SFR values for this particular galaxy are not unex- 
pected. 

The sub-mm (e.g. 350 GHz) emission from galaxies serves 
as another estimator of SFR, and it is related to the radio emis- 
sion at 1 .43 GHz via a reds hift-dependent spectral index, /? 35 P 
(|Carilli & Yun 1999|, 2000; |Dunne, Clements & Eales 2000|). 



Using the calibration of Dunne, Clements & Eales (2000) we 
find a value of /3 35 4 ° 0.6 ± 0.1 at z « 1, which gives F„(350) 
1-9^'g mJy. Observations with the Sub-millimeter Common 
User Bolometer Array (SCUBA) camera on the James Clark 
Maxwell Telescope (JCMT) 12.4 days after the burst provided 
a 2a upper limit of 3.2 mJy on t he combined emissi on from the 
afterglow and host at 350 GHz ( Smith et al. 1999 ), consistent 
with the predictions from the radio-sub-mm relation. 

To conclude, we use the derived SFR to calculate the ex- 
pected far-IR (FIR) emission from the host of GRB 980703. 
The luminosity of the FIR radiation can be derived from the em- 
pirical relation suggested by Helou, Soifer & Rowan-Robinson 
(1985), 

q = -12.6+log(F FIR /F 1A ) » 2.3 (2) 
1O I2 L for the host of GRB 980703; 



which evaluates to Lf /« : 
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erg sec 1 cm 2 



here Ffw is the total flux in the range 40-120 in units of erg 
and F\ .4 is the flux density at 1.4 GHz in units of 
Hz -1 . Hopkins et al. (2001) provide the calibra- 
tion L Fm = 5.81 x 10 9 x SFR L Q , which gives L FIR w 3 x 10 12 
L©, in good agreement with the previous value. These values 
of the FIR luminosity pl ace the host galaxy of GR B 980703 in 
the category of ULIRG ( [Sanders & Mirabel 199^ ). 



burst, give values of 1-2 mags from optical observations, and 
som ewhat higher values from X-ray observations (Bloo m et al. 



5. OFFSET MEASUREMENTS AND SOURCE SIZE 

Figure [| shows the projected angular offset between the host 
galaxy and afterglow of GRB 980703, for each individual de- 
tection, and the combined value from all observing runs (insert 
in Figure Q). Positions are determined from Gaussian fits, and 
the host-GRB offset is calculated with respect to a Very Long 
Baseline Array (VLBA) position that was measured to 0.0007 
arcsec accuracy in each coordinate on 1998 August 2 at 8.42 
GHz (Bergeretal. 2001). 

We find an average offset from all measurements of -0.032 ± 
0.015 arcsec in RA and 0.024 ±0.015 arcsec in declination. 
The uncertainty in the position of the source is given by 
(^offset ~ (6 S y„beam/2)/(S/N), where Qsynbeam ~ A/fi„ MA - is the 
half-power synthesized beam-width, A is the observing wave- 
length, B max is the length of the maximum baseline, and S/N is 
the signal-to-noise ratio of the flux measurement. 

The optical measurements of Bloom, Kulkarni & Djorgov- 
ski (2000) for the host of GRB 980703, give an angular offset 
of -0.054 ± 0.055 in RA and 0.098 ± 0.065 in declination (see 
insert in Figure ||). They conclude that GRB 980703 was not 
significantly offset from the center of its host galaxy, in agree- 
ment with the more accurate offset measurements in the radio. 

In addition to accurate measurements of the offset, the ra- 
dio observations allow us to place meaningful limits on the size 
of the radio-emitting region (i.e. the size of the star-forming 
region). We find that in our highest resolution images the 
source is unresolved, and therefore, based on the synthesized 
beam size we can derive an upper limit on the physical size 
of the source. For our adopted cosmological parameters (sec- 
tion we find that the angular diameter distance to the source 
is (1a ~ 5.4 x 10 27 cm. The full synthesized beam-width at 8.46 
GHz is 9hpbw ~ 0.27 arcsec, which gives an upper limit of 
Dmd = AaGhpew < 2.3 kpc on the diameter of the source. 

Holland et al. (2001) used a R l/n profile to fit the optical 
emission from the host and found that the best fit gives a half- 
light radius of 0. 1 3 arcsec, which corresponds to an exponential 
disk with a scale-diameter of 0.44 arcsec. Thus, the physical 
size of the galaxy is D opt « 3.7 kpc, 60% larger than the up- 
per limit from our radio measurements. However, Holland et 
al. (2001) claim that the center of the galaxy is 0.2 mag bluer 
than the outer regions of the host. If so, star formation must be 
mainly taking place within the inner parts of the galaxy. Since 
the radio emission directly traces current star formation, we ex- 
pect the radio emission to be more centrally concentrated than 
the optical emission. Thus, as expected, the radio size of the 
galaxy is smaller than the optical size. 

Most likely the GRB is located within the nuclear starburst 
given the small offset of the GRB from the centroid of the 
galaxy. If so it raises the question of why the afterglow was not 
completely extinguished by dust. In fact, in order to reconcile 
the optical and radio derived SFRs we require a rest-frame ex- 
tinction of Ay ~ 4.5 mag. Observations of the afterglow, which 
provide an estimate of extinction along the line-of-sight to the 



1998; |Castro-Tirado et al. 1999| ; [Vreeswijk et al. 1999[ ). Thus, 
the extinction in the nuclear star-forming region is higher than 
the average over the whole galaxy, and the correction to the ob- 
served optical SFR is almost sufficient to reconcile it with the 
value of 500 M Q /yr derived from the radio. 

The relatively small source size also agrees well with the 
classification of the host of GRB 980703 as a ULIRG exhibiting 
a starburst. Kennicutt (1998, and references therein) showed 
that star formation with a rate ;> 20 M©/yr invariably takes 
place in circumnuclear regions of size 0.2-2 kpc, in the form 
of nuclear starburst. As a result, we expect that ULIRGs will 
have such size scales when traced by star formation, and the 
source size we measured for the host of GRB 980703 indicates 
that it is probably undergoing a nuclear starburst. 

Finally, from the source size and offset measurement we 
conclude that GRB 980703 took place in the region of maxi- 
mum star formation, providing further indirect evidence linking 
GRBs to massive stars 

6. CONCLUSIONS 

Late-time observations of GRB 980703 reveal a steady com- 
ponent, with a flux density = 68.0 ±6.6 /iJy and a spec- 
tral index (3 = 0.32 ±0.12. The spectral and temporal charac- 
teristics of this emission indicate that it does not arise from the 
afterglow itself, but rather it is the result of star formation in the 
host galaxy, with SFR « 500 M Q /yr. This leads to the interpre- 
tation that this host galaxy is a ULIRG undergoing a starburst. 
In addition, the star formation is concentrated within the inner 
two kpc of the host, and the progenitor of GRB 980703 was po- 
sitioned within this region of star formation. This conclusion 
lends additional support for the collapsar model. 

If GRBs really come from massive stars, then they can be 
used to trace the star for mation history of the universe (e.g. 
Blain & Natarajan 20"00| ). In addition, GRBs and their af- 
t erglows are potentially detectable out to very high redshifts 
(Lamb & Reichart 200C). These propositions, taken together 



with the dust-penetrating power of their 7-ray emission, make 
GRBs a unique tool for the study of galaxies and the IGM over a 
wide redshift range. In particular, radio and sub-mm/FIR obser- 
vations of a GRB-selected galaxy sample will be extremely use- 
ful for the study of the obscured star formation fraction, and the 
properties of starbursts at high redshifts. Moreover, a compar- 
ison of the global star formation history as derived from these 
long-wavelength host studies, with the redshift distribution of 
GRBs, will provide valuable insight as to how well GRBs trace 
the formation rate of massive stars; we expect that if GRBs trace 
only a particular channel of star formation, the two distributions 
will not agree. 

Therefore, it is imperative to study the hosts of GRBs in the 
radio and sub-mm/FIR. Future observatories such as the Space 
Infrared Telescope Facility (SIRTF; to be launched in July 
2002), the Expanded VLA (EVLA), and the Square-Kilometer 
Array (SKA) will allow detailed studies of these hosts. In the 
FIR, SIRTF will have the ability to detect sources down to a 
few mJy, allowing the detection of galaxies with SFR compa- 
rable to that in the host of GRB 980703 out to z ~ 10; alter- 
natively, we will be able to detect hosts with SFR as low as 
a few M Q /yr at z ~ 1 . The EVLA and SKA will greatly im- 
prove the detectability of host galaxies in the radio, and will 
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also allow much higher angular resolution studies of compact 
star-forming regions. With a factor ten increase in resolution 
and a factor five increase in sensitivity over the current VLA, 
we will be able to probe scales of approximately 5 mas with the 
EVLA; for a galaxy at z « 1 this translates to a physical scale 
of 150 pc. In addition, EVLA will detect galaxies with a total 
SFR as low as 50 M Q /yr at z ~ 1. The SKA, with a similar 
resolution, but a much larger collecting area, will extend this 



capability to even lower SFR and smaller star-forming regions. 

Thus, as more host galaxies are detected and studied in detail 
in the radio and sub-mm/FIR, we will be able to address a large 
number of issues pertaining not only to the bursts themselves, 
but also to the characteristics of galaxies at high redshifts. 

We acknowledge support by NSF and NASA grants. 
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Table 1 

Radio Observations of GRB 980703 



Epoch 


At 


Array Configuration 


^on-source 


^0 


S±CT 


(UT) 


(days) 




(hrs) 


(GHz) 


(MJy) 


1999 June 15.36—26.29 


352.65 


A 


6.5 


1.43 


81±18 


1999 July 10.53—28.28 


381.22 


A 


13.4 


1.43 


57±10 


1999 August 19.40— September 21.24 


428.64 


A 


11.3 


8.46 


37±8 


1999 November 24.06 


508.88 


B 


1.7 


1.43 


57±20 


2000 March 5.70 


610.52 


BnC 


2.8 


8.46 


57±14 


2000 October 7.30— November 19.08 


847.01 


A 


6.5 


1.43 


76=1=11 


2000 December 2.19 — 4.97 


882.60 


A 


4.9 


4.86 


35±11 


2000 December 21.15— 2001 January 4.98 


908.38 


A 


7.0 


8.46 


41±8 


2001 February 2.00^1.93 


945.29 


AnB 


1.9 


1.43 


113±21 


2001 February 8.08 


949.90 


AnB 


1.2 


4.86 


57±28 


2001 March 2.98—9.00 


975.81 


B 


4.5 


4.86 


43±15 


2001 March 22.96— April 8.56 


1001.08 


B 


2.6 


1.43 


83±30 



NOTE. — The columns are (left to right), (1) UT date of the start of each observation or range of dates for observations which were added over several days, (2) time 
elapsed since the 7-ray burst, (3) array configuration, (4) total on-source observing time, (5) observing frequency, and (6) peak flux density at the best fit position of 
the radio transient, with the error given as the root mean square noise on the image. 
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Fig. 1. — Radio lightcurves at 1.43, 4.86, and 8.46 GHz. The thin solid lines are the combined afterglow and host galaxy emission, 
the dotted lines indicate the afterglow emission, and the thick solid lines are the weighted-average fluxes of the host galaxy, with the 
thickness indicating the uncertainty in the flux. Only measurements at t ;> 350 days after the burst were used to calculate the host 
flux. The fits are based on broadband fitting (Berger et al. 2001). The data clearly indicate that there is a constant component in the 
observed emission, interpreted as the host galaxy. For details of the data reduction and analysis see section [| 
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Fig. 2. — Fluctuations of individual measurements around the weighted average presented by the wide strips in figure [j] We note 
that there are no fluctuations above l.ler at 8.46 GHz, 0.8er at 4.86 GHz, and 1.7er at 1.43 GHz, indicating that the flux in each band 
is consistent with a constant. In fact, if we assume that the source has some variability (over that due to measurement errors), then 
the variable flux is less than ±7 /iJy, with a 40% probability that there are no fluctuations at all (see insert). This conclusion supports 
the hypothesis that the radio emission is not due to an AGN. 
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Fig. 3 — SFR for radio galaxies in the HDF, SSA13 and V15 (Haarsma et al. 2000) in the redshift range 0.85-1.15. The dashed 
line is an average over this redshift range. The host galaxy of GRB 980703 is marked by a solid square. The errorbars on the survey 
sources indicate a 10 niy measurement error. The host of GRB 980703 appears to be an average star-forming galaxy relative to the 
sample. 
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Fig. 4. — Offset measurements for all epochs in which the host galaxy emission dominates. The plot shows the offset in RA and 6 
of the VLBA position of the burst (see section |[) relative to the host center, (A#ra, A9s)=(0,0). The most accurate measurements 
are at 8.46 GHz in the VLA A configuration. In this mode we achieved an rms positional error of 0.02 arcsec. The insert shows 
the weighted average offset in both RA and S (small cross). The larger cross is the offset measurement from Bloom, Kulkarni, and 
Djorgovski (2001). The solid circle designates the projected maximum source size from the radio observations in the A configuration 
at 8.46 GHz, and the dashed circle is the optical size from Holland et al. (2001). Clearly the formation of massive stars is concentrated 
in the central region of the host. The small offset of the burst from the host center indicates that GRB 980703 occurred in the region 
of maximum star formation, which points to a link between GRBs and massive stars. 



